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Los Alamos Meson Physics Facility by using from fj, + decay 
at rest. The i> e are detected via the reaction v e p — > e + n, 
correlated with a 7 from np — * d"f (2.2 MeV). The use of 
tight cuts to identify e + events with correlated 7 rays yields 
22 events with e + energy between 36 and 60 MeV and only 
4.6 ± 0.6 background events. A fit to the e + events between 
20 and 60 MeV yields a total excess of 51.8±i|; % ± 8.0 events. 
If attributed to Dp — » oscillations, this corresponds to an 
oscillation probability of (O.Sl^J ± 0.05)%. 
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We present the results from a search for neutrino os- 
cillations using the Liquid Scintillator Neutrino Detector 
(LSND) apparatus described in reference |Q. The exis- 
tence of neutrino oscillations would imply that neutrinos 
have mass and that there is mixing among the different 
flavors of neutrinos. Candidate events in a search for the 
transformation — > v e from neutrino oscillations with 
the LSND detector have previously been reported || for 
data taken in 1993 and 1994. Data taken in 1995 have 
been included in this paper, and the analysis has been 
made more efficient. 

Protons are accelerated by the LAMPF linac to 800 
MeV kinetic energy and pass through a series of targets, 
culminating with the A6 beam stop. The primary neu- 
trino flux comes from ir + produced in a 30-cm-long water 
target in the A6 beam stop jjj. The total charge deliv- 
ered to the beam stop while the detector recorded data 



was 1787 C in 1993, 5904 C in 1994, and 7081 C in 1995. 

Most of the 7T + come to rest and decay through the 
sequence n + — ► [i^u^, followed by /i + — » e + ^ e j> M , supply- 
ing Dp with a maximum energy of 52.8 MeV. The energy 
dependence of the flux from decay at rest (DAR) is 
very well known, and the absolute value is known to 7% 
plpj. The open space around the target is short com- 
pared to the pion decay length, so only 3% of the 7r + 
decay in flight (DIF). A much smaller fraction (approxi- 
mately 0.001%) of the muons DIF, due to the difference 
in lifetimes and that a ir + must first DIF. The total P M 
flux averaged over the detector volume, including con- 
tributions from upstream targets and all elements of the 
beam stop, was 7.6 x 10 _10 P A1 /cm 2 /proton. 

A v e component in the beam comes from the symmet- 
rical decay chain starting with a 7r~. This background 
is suppressed by three factors in this experiment. First, 
7r + production is about eight times the 7r~ production 
in the beam stop. Second, 95% of 7r~ will come to rest 
and are absorbed before decay in the beam stop. Third, 
88% of fi~ from 7r~ DIF are captured from atomic orbit, 
a process which does not give a v e . Thus, the relative 
yield, compared to the positive channel, is estimated to 
be ~ (1/8) x 0.05 x 0.12 = 7.5 x 10~ 4 . A detailed Monte 
Carlo simulation HI, gives a value of 7.8 x 10~ 4 for the 
flux ratio of v e to v^. 

The detector is a tank filled with 167 metric tons of 
dilute liquid scintillator, located about 30 m from the 
neutrino source, and surrounded on all sides except the 
bottom by a liquid scintillator veto shield. The dilute 
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mixture allows the detection in photomultiplier tubes 
(PMTs) of both Cerenkov light and isotropic scintilla- 
tion light, so that reconstruction provides robust particle 
identification (P1D) for along with the position 
and the direction of the event. PID is based on the qual- 
ity of the position and Cerenkov angle fits, and on the 
relative amount of early light Q|. The detector needs 
to distinguish between events induced by T> e (oscillation 
candidates) from the events produced by the v e . LSND 
detects v e via u e p — > e + n, a process with a well-known 
cross section [0 , followed by the neutron-capture reaction 
np — > d 7(2.2 MeV). Thus the oscillation event signature 
consists of an "electron" signal, followed by a 2.2 MeV 
photon correlated with the electron signal in both posi- 
tion and time. Detection of DAR v e in LSND is domi- 
nated by charged current reactions on 12 C, but an elec- 
tron from v e 12 C — > e~ 12 N has energy E e < 36 MeV be- 
cause of the mass difference of 12 C and the lowest lying 
12 N state. Moreover, the DAR production of a corre- 
lated photon from v e 12 C — > e~n 11 N can occur only for 
E e < 20 MeV because of the threshold for free neutron 
production. 

Cosmic rays are suppressed at the trigger level by use of 
the veto shield and by rejecting events with any evidence 
for a muon in the previous 15.2 fis 0j. Even so, the 
trigger rate is dominated by this background, with actual 
^-induced events contributing less than ~ 10~ 5 of all 
triggers. Because the data acquisition and triggering 
do not depend on whether the beam is on or off, the 
beam-on to beam-off duty ratio could be measured from 
triggered events; it averaged 0.070 over the three years of 
data. The beam-unrelated background in any beam-on 
sample is thus well measured from the much larger beam- 
off sample, and can be subtracted. The cuts used to 
select e + candidates are designed to discriminate heavily 
against this background, so that the statistical error from 
this subtraction can be kept small relative to the beam- 
dependent signal. 

Separation of correlated neutron-capture photons from 
accidental signals is achieved using an approximate like- 
lihood ratio, R MM, f° r the correlated and accidental 
hypotheses. R is defined using distributions || of the 
number of hit PMTs for the reconstructed 7 and of the 
time and distance between the primary event and that 
7. For purposes of fitting, the R distribution for acciden- 
tal photons is taken from laser calibration events. That 
for correlated photons is taken from cosmic ray neutron 
events either directly or as modified for the lower-energy 
neutrons of interest by using a Monte Carlo simulation of 
the distance distribution, with fit results averaged over 
the two cases. 

We present analysis of the full 1993+1 994+1995 data 
sample for two sets of positron selection cuts. Selection A 
corresponds to the criteria used in our previous paper on 
the 1993 and 1994 data ||. Selection B uses new insight 
into the nature of the beam-off backgrounds to further 



reduce these backgrounds while relaxing other criteria to 
increase the signal efficiency by about 40%. The crite- 
ria were chosen, and efficiencies determined, using several 
control samples taken as part of the data stream. A sam- 
ple of "Michel" electrons from the decays of stopping cos- 
mic ray muons is used to characterize energy calibration, 
resolution and PID. Cosmic ray neutrons stopping in the 
detector are used for the 2.2 MeV 7 properties and as 
a "non-electron" control sample for electron PID. Other 
neutrino induced interactions in the detector including 
12 C — > [i~X 0, and v e 12 C — > e~X are also used to 
check efficiencies and backgrounds. Random triggers in 
association with tank calibration are used to determine 
veto efficiencies, readout deadtime, and the distribution 
of R for accidentally coincident 7. 

The primary particle in a v e event candidate is required 
to have a PID consistent with a positron. The selection A 
criteria for PID were previously described ||], giving an 
efficiency for positrons in the 36 < E e < 60 MeV energy 
range of 0.77 ±0.02. Selection B loosens the PID criteria 
to increase the positron PID efficiency to 0.84 ± 0.02. 

Selection A removed all events with the time to the 
previous triggered event, At p < 40 (is to eliminate Michel 
electrons from muon decay. Selection B required At p 
greater than 20 /is, and no activities between 20 (is and 

34 /is before the event trigger time with more than 50 
PMT hits or reconstructed within 2m from the positron 
position. The selection A and B efficiencies are 0.50±0.02 
and 0.68 ±0.02, respectively. The time to any subsequent 
triggered event, At a , is required to be > 8/iS (~ 4 muon 
lifetimes) to remove events which are misidentified muons 
which decay (0.99 ± 0.01 efficiency). The reconstructed 
positron location was required to be a distance D > 35 
cm from the surface tangent to the faces of the PMTs 
(0.85 ± 0.05 efficiency). This assures that the positron 
is in a region of the tank in which the energy and PID 
responses vary smoothly and are well understood. The 

35 cm cut also avoids the region of the tank with the 
highest cosmic ray background. 

To suppress cosmic ray neutrons, the number of asso- 
ciated 7 with R > 1.5 is required to be no more than 
2 for selection A (0.99 ± 0.01 efficiency) and no more 
than 1 for selection B (0.94 ± 0.01 efficiency). Cosmic 
ray neutrons that enter the detector often produce one 
or more additional neutrons, while recoil neutrons from 
the v e p — > e + n reaction are too low in energy to knock 
out additional neutrons. The number of veto shield hits 
associated with the events is no more than 1 for selection 
A (0.84 ±0.02 efficiency) and no more than 3 for selection 
B (0.98 ±0.01 efficiency). 

Beam-off data surviving these cuts differ from the ex- 
pected neutrino interaction signal in two respects. One is 
the distribution of angles between the e + direction and 
its position vector relative to the tank center - back- 
ground events tend to head inwards. The other is in the 
distribution of veto hits - cosmic ray events tend to have 
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TABLE I. The number of signal and background events in 
the 36 < E e < 60 MeV energy range. Excess/Efficiency is the 
excess number of events divided by the total efficiency. The 
beam-off background has been scaled to the beam-on time. 
B' is a restrictive geometry test. 



Selection 


Signal 


Beam-Off 


v Bkgd. 


Excess 


Excess 
Efficiency 


A R > 


221 


133.6 ± 3.1 


53.5 ± 6.8 


33.9 ± 16.6 


130 ± 64 


A R > 30 


13 


2.8 ± 0.4 


1.5 ± 0.3 


8.7 ± 3.6 


146 ± 61 


B R > 


300 


160.5 ± 3.4 


76.2 ± 9.7 


63.3 ± 20.1 


171 ± 54 


B _R. > 30 


22 


2.5 ± 0.4 


2.1 ± 0.4 


17.4 ± 4.7 


205 ± 54 


B' R > 


99 


33.5 ± 1.5 


34.3 ± 4.4 


31.2 ± 11.0 


187 ± 66 


B' R. > 30 


6 


0.8 ± 0.2 


0.9 ± 0.2 


4.3 ± 2.5 


110 ± 63 



more of them. These two distributions are used in a way 
analogous to the R parameter discussed earlier in defin- 
ing a likelihood ratio, S ||, which is used as the final 
positron selection criterion. For selection B, but not A, 
we require S > 0.5, a cut that loses 13% of the expected 
neutrino signal while eliminating 33% of the beam-off 
background. Including a 0.97 ±0.01 data acquisition effi- 
ciency gives overall efficiencies of 0.26 ±0.02 for selection 
A, and 0.37 ± 0.03 for selection B. 

The backgrounds to v e p — > e + n followed by n capture 
fall into three general classes: beam-off events (cosmic 
ray induced), beam- related events with correlated neu- 
trons, and beam-related events with an accidental 7. As 
outlined above, the cosmic ray background to beam-on 
events is 0.07 times the number of beam-off events which 
pass the same criteria. 

The major sources of beam-induced backgrounds are 
from yT DAR, discussed above, and from 7r + DIF in 
the beam stop. The latter results in a background from 
p interactions where the final /i + is missed, and its 
Michel decay positron is mistaken for a primary v e p 
event. These backgrounds are estimated using the 
detector Monte Carlo simulation The backgrounds 

with accidental 7 overlap are greatly reduced by selection 
on the R parameter. Details of all backgrounds consid- 
ered are presented in Rcf . |j| . 

Table I lists the number of signal, beam-off background 
and neutrino-background events for the two selections 
with 36 < E e < 60 MeV - to avoid large accidental-7 
backgrounds. The likelihood ratio, R, is used to deter- 
mine whether a candidate 2.2 MeV 7 is correlated with 
an electron or from an accidental coincidence. Requir- 
ing R > 30 (correlated-7 efficiency = 0.23) we observe 22 
events beam-on and 36 x 0.07 = 2.5 events beam-off. The 
estimated beam-related background consists of 1.72±0.41 
events with correlated neutrons, and 0.41 ±0.06 without. 
The probability that the beam-on events are entirely due 
to a statistical fluctuation of the 4.6 ±0.6 event expected 
total background is 4.1 x 10~ 8 . Figure [j](a) shows the 
energy distribution of all primary electrons which pass 
selection B with associated R > 0. Figure 0(b) shows the 
electron energy distribution for selection B with R > 30. 

Kolmogorov tests have been done to check for unex- 
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FIG. 1. The energy distribution for events which pass se- 
lection B with (a) R > and (b) R > 30. Shown in the fig- 
ure are the beam-excess data, estimated neutrino background 
(dashed), and expected distribution for neutrino oscillations 
at large Am 2 plus estimated neutrino background (solid). 



pected concentrations of events in position (e.g., in re- 
gions of high cosmic ray or 7 backgrounds), energy or 
time (year). No consistency check yields a probability 
so low as to demonstrate a serious inconsistency ||. A 
restrictive geometric cut, removing the 55% of the selec- 
tion B acceptance with highest cosmic ray rates [|| , also 
demonstrates no inconsistency; its results are labelled B' 
in Tabic | 




FIG. 2. The R distribution, beam on minus beam-off ex- 
cess, for events that satisfy selection B and that have energies 
in the range 20 < E e < 60 MeV. The solid curve is the best fit 
to the data, the dashed curve is the uncorrelated 7 component 
of the fit, and the dotted curve is the correlated 7 component. 

To determine the oscillation probability we fit the over- 
all R distribution, for events satisfying selection B, in the 
full energy range 20 < E e < 60 MeV. The larger energy 
range is used in this and the following fit to utilize the 
maximum amount of data and is made possible by our in- 
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creased understanding of the background processes. The 
1763 beam-on and 11981 beam-off events were fit by a x 2 
method which took spatial variations in accidental pho- 
ton rates into account by averaging the appropriate R 
distributions at the positions of each positron. The result 
of the fit is shown in Fig. [| It yielded 64.3l}g;y beam- 
related events with a correlated 7, and 860^g beam- 
related events without a correlated 7. The latter is con- 
sistent with a calculated background estimate of 795±134 
such events. Subtracting the estimated neutrino back- 
ground with a correlated 7 (12.5 ±2.9 events) results in a 
net excess of 51.8 events, corresponding to an oscilla- 
tion probability of (0.3ll ; 10 ± 0.05)%, where the second 
error is systematic. A likelihood fit which uses individual 
local accidental-7 R distributions for each positron gave 
a consistent result of (0.27t£;:[ 2 ± 0.04)%. 
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FIG. 3. Plot of the LSND Am 2 vs sin 2 20 favored regions. 
The shaded regions are the 90% or 99% likelihood regions as 
defined in the text, not confidence regions. Also shown are 
90% C.L. limits from KARMEN at ISIS (dashed curve), E776 
at BNL (dotted curve), and the Bugey reactor experiment 
(dot-dashed curve). 

For simplicity we present the results in the two- 
generation formalism, in which the mixing probability 
is written as P — sin 2 20sin 2 (1.27ATO 2 L/^), where 6 is 
the mixing angle, Am 2 is the difference of the squares of 
the two mass eigenstates in eV 2 , L is the distance from 
neutrino production in meters, and E v is the neutrino 
energy in MeV. An overall likelihood fit has been made 
to determine favored regions in the Am 2 versus sin 2 Id 
parameter space for two-neutrino mixing. The fit was 
made to distributions in the observed event energy, the 
neutron likelihood ratio R, the reconstructed direction of 
the electron relative to the neutrino beam direction, and 
the distance of the primary event from the beam stop 
neutrino source. The beam-related and cosmic ray back- 



grounds were added to the expected neutrino oscillation 
signal, and a likelihood was calculated for a range of Am 2 
versus sin 2 29 values. Figure || shows regions which are 
within 2.3 and 4.5 log-likelihood units of the maximum, 
called 90% or 99% likelihood regions. The regions have 
been enlarged to account for systematic effects by varying 
the inputs to the fit to reflect uncertainty in backgrounds, 
neutrino fluxes and the R distribution shape. Figure 
also shows the 90% C.L. limits from KARMEN at ISIS (f 
(dashed curve), E776 at BNL (dotted curve) |§, and the 
Bugey reactor experiment Q (dot-dashed curve). 

This paper reports the observation of 22 electron 
events in the 36 < E e < 60 MeV energy range that 
are correlated in time and space with a low-energy 7 
with R > 30, and the total estimated background from 
conventional processes is 4.6 ± 0.6 events. The prob- 
ability that this excess is due to a statistical fluctu- 
ation is 4.1 x 10~ 8 . A fit to the full energy range 
20 < E e < 60 MeV gives an oscillation probability of 
(0.31±£ u ± 0.05)%. These results may be interpreted 
as evidence for — > D e oscillations within the favored 
range of Fig. S. 
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